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1 .0 INTRODUCTION 


lha fifth monthly prograss raport includas corractions and additions 
to tha praviously subaittad reports. Iha addition of tha SRB propallant 
thicknass as a stata variabla is incl\ided with the associated partial 
derivatives . 

Curing this reporting period j praliminaury results of tha estimation 
program checkout was presented to NASA technical personnel. 


i 


Ihe system described by equation <1 ) is observed at discrete times > t^^i 
with not all states being directly measured. Some measurements are non- 
linear functions of the elements of the state vector x(t). In general 
the measurement process is described as 

( 2 ) 

%d»ere 

^ * m-dimensional observation vector 

h^ 3 functional representation of the measurements in terms of 
the states 

3 m-dimensional I zero-mean > lAile Gaussian noise sequence with 
covariance 

E[v. vT] = R. 6. . 

-1 “3 1 1*3 

Examples of the elements of the observation vector ^ include radar 
measurements of range, azimuth, and elevation from the radar site to the 
vehicle. 

It is assumed that the system process noise vector w(t) and the 

measurement noise vector ^ are uncorrelated. Also, the system state 

initial condition vector x is not correlated with either of these two 

-o 

noise vectors. Therefore 

E[w(t) v^] 3 0, E[w(t) x"^] 3 0, E[x yj] s 0 
— -o -o He 

where the superscript T denotes transpose. Fox later reference, the 
following matricies are defined 


2.0 FILTERING AND SMOOTHING ALGORITHM 


Die Space Shuttle Parameter Estimation Program utilizes optimal 
estimation techniques to provide estimates of the propulsion system 
peurameters . Ihe technique selected is the extended Kalman filter end 
the modified Bryson -Frazier onoother. By modeling the propulsion system 
parameters as time correlated remdom variables , improved estimates of 
t^-^se parameters are obtained and are properly time phased by removing 
the filter induced lag by using the combined f liter /smoother . The 
smoother also provides improved estimates of the initial state estimates 
Hie system I in state-space notation, is modeled as the continuous 
dynamical system disturbed by additive Gaussian vAiite noise 

X « , t ) -»■ G(t) w(t) + ^(t) , x(o) ■ Xq ) 

where 

x > n -dimensional state vector 

x a Gaussian initial condition vector with covauriance P 
-o o 

w(t) a p-dimensional white, zeri-mean vAiite Gaussian noise with 
covariance 

E[w(t) w^(t)] a Q(t) 6(t - t) 

^(t) a n-d-ijnensional control vector. 

The elements of the vector x(t) represent vehicle position, velocities 
attitudes, angular rates, aerodynaimic and propulsion parameters, measure 
ment biases, etc. Elements of u(t) include known control inputs such 


as SSME power level commands . 



# ’ w 


af(x(t), t) 

F^x(t), t) - ■. (3) 

and 

ah (x (t ) ) 

— k 

hese matricies are linearizations of the dynamics and measure.nent models 
respectively, evaluated about either a nominal or reference value of the 
state, or about the state estimate. 

2.1 Extended Kalman Filter Algorithm 

Ihe extended Kalman filter al<jorithm is in essence a conventional 
linear Kalman filter algorith*" applied to a mathematical model resulting 
from the linearization of the system model equation (1), and measurement 
process, equation (2), about a current state estimate. Hie filter yields 
optimal estimates if the linearization is accurate, i.e., the state esti- 
mate closely approximates the tr»ie state. The derivation of the algorithm 
can be found in reference [1 ] . 

The algorithm proceeds as follows. After initialization of the state 
estimate and covariance, the state estimate and covsuriance sure propagated 
forward in time until a measurement update is available, by 

X = ^(x(t) , t ) , ^k 1 ^ ^ “ ^k ^ 

and 

P(t) « F(x(t), t) P(t) + P(t) F(x(t), t)"^ + G{t) Q(t) GCt)"^ 


( 6 ) 


At the measurement time, the state estimate and covariance are updated 


by 


and 

P^( + ) - (I - P^(-) (8) 

where the (-) and (+) represent the appropriate values just before 
and just cifter the update. The updated values are used to reinitialize 
the time propagation equations (3) and (4) for integrating up to the next 
measurement time. The Kalman gain matrix is computed as 

This algorithm is repeated until the last time point, t^^, is processed. 

For later use in the smoother algorithm, various combinations of the state 
estimates (x), measurements (£), linearized dynamics matrix (F) and 
measurement matrix (H), measurement noise covariance (R) and estimation 
error covariance matrix (P) must be stored for each time instant to be 
processed by the smoother algorithm. 

2.2 Modified Bryson -Frazier anoother Algorithm 

The operation of the smoother algorithm is similar to the filter 
algorithm except in reverse time. The derivation of this smoother algorithm 
is found in reference [2]. This fixed interval smoothing algorithm pro- 
vides optimal estimates given all the measurements in comparison to the 
filtering algorithm providing optimal estimates given the previous 


measurements processed. ‘Rierefore the smoother provides Improved estimates 
in addition to removing the time lag induced by the filter algorithm. 

Ihe smoothing algorithm ^joint variables » ^ and A are "initialized" 
at the final time processed by the filter > T, 

X(T-) . -hJ(H hI . (10) 

N,T 

and 

- «n‘“n '’n "n * %>■' "n V , 

N,T 

If T is not an observation time. ^ and A are zero. The ^ joint 
variables are propagated in reverse time to the next previous measurement 
time by 

X - - F(x(t), t) X , 1 t < tj^^^ (12) 

A - - F(x(t), t)'^ A - A F(x(t), t) (13) 

At the time of an available measurement) t^^, the adjoint variables 
are updated by 

X(-) - X(+) - H^(H^ Pfc 

+ (t^ -t- \) kJ X( + )) (14) 

2Uld 

. (I - \ "k * V' \ <' = ' 


A(-) 


The smoother state estimate and error covariance are obtained using the 
filter estimate and covariance and the adjoint variables by 

x*(t) - kit) - P(t) X(t) (16) 

and 

P*(t) « P(t) - P(t) A(t) P(t). (17) 

Due to the potential number of time points to be processed, smoother 
estimates may only be ccmtputed at the discrete measurement times. For 
this approach the propagation equations (10) and (11) are replaced by 

(18) 

and 

+ ) A. , (-) ♦,, (19) 

k k k+1 k 

T 

where 4, is the state transition matrix formed with the linearized 
k 

dynamics matrix F to propagate the adjoint variable from time 
to time t^. The algorithm continues in reverse time until the initial 
time is reached. 

2.3 Iterations with the Filter /awoother Algorithm 

The performance of the filter/smoother algorithm is a direct result 
of the accuracy of the linearization. Repeated operations of the algorithms 
with adjustments in initial state estimates and covariance in each cycle 


can yield improved estimates. This technique is known as global iterated 


filtering as defined in reference (31. Each cycle of operating the 
algoritluns would yield increasing improvements in the state estimates. 

This feature of the algorithm operation is of special interest to 
the propulsion petrameter estimation problem using the NASA predictive 
models. Initial I or nominal « values of the parameters of interest can 
be used to obtain the necessary partial derivatives indicated earlier. 
From operating the algorithm improved estimates of those parameters are 
obtained. Using these improved estimates* mors accurate partial deriva- 
tives ore obtained for use in the algorithms. This process is continued 
until there is in essence no change in the partial derivatives or quality 
of the state estimates. If the linearization is accurate, the measure- 
ment residual should be a vdiite noise process with known covariance. 


3.0 FILTER/SMOOTOER ALGORITHM SYSTEM AND MEASURgMENT MODEL 


Iha utafulnesB of the f iltar/smoother algorithm la to provide esti- 
mates of the system states fr<Mn the observed motion and dynamics while 
the syst«n is driven by known and unknown elwnents. 'n\ese unknown elements 
are elements of the system state vector to be estimatec.. Ihe evolution of 
motion resulting from these known and unknown elements is assumed to be 
suitably represented for this study by a six degree -of -freedom (6 DOF) rigid 
body equations of motion. These equations are presented and discussed in 
section 3.1. 

To implement these equations into the filter /smoother algorithm 
presented in section 2.0, a linearization nf the system state and measure- 
ment models is required. These linearized equations are presented in 
section 3.2. 


3 .1 Equations of Motion and Measurement Equations 


3.1.1 Rigid Body Equations of hotion 

The rate of change of vehicle velocity in body coordinates, 

as a result of external forces acting on the vehicle is described 
2 


(B) 


pAv 

2ra 


(B) (B) 

m B^I (I) , (I), (B) if 

c_+Cg (r )-(uxv + + — — — 

-f *■ — — — mm 


(B) 

V , 

by 


( 20 ) 


where 

p ■ atmospheric density 

A - aerodynsnic coefficient referenced area 

V ■ magnitude of vehicle velocity relative to the surrounding air 


mass 


m ■ Vtthicle mas* 


^ ■ aarodynamlc force coefficient vector 

^(I) _ gravity vector in inertial coordinate* 

^ • angular rotation of the body relative to the inertial frame 
(Q ) 

^ “ resultant thrust force vector in body coordinates 

(B) 

fj, ■ resultant plvkae force vector in body coordinates 

The rate of change of vehicle position in inertial coordinates, 
is then obtained by 


*(I) IB (B) 
r *» C V 


(21 ) 


where C is the transformation matrix from body coordinates to 

I B 

inertial coordinates. The elenents of the C tr an c formation matrix 

are obtained from the resulting Euie,. wJi^les defined by 
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1 sinftanO ccsgtani 


P 

• 

e 

- 

0 cosy -sln9 


q 



0 sinfsecO cos^secB 


r 

» - 




» - 


v^ere 9 , 6 , and are roll, pitch and yaw attitudes respectively. 

The roll, pitch and yaw rates of the body relative to inertial coordinates 
are p, q, and r respectively. Finally, the rate of change of the body 
rates relative to inertiil is given by 
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(B), 

r ) X c. 
-eg -f 


- Cl) X (^o>) + T„ + T„ J 
— — -T P 


where 


I = vehicle moments of inertia matrix 

= aerodynamic coefficient referenced Ic.ngth and mcxnent 

coefficient vector 
( B ) 

r^^ = vehicle center -of -gravity vector in body coordinates 

r^ = aerodynamic coefficient reference position in body coordinates 
( B ) 

^ = resultant thrust torque vector in body coordinates 

( B ) 

= resultant plane torque vector in body coordinates 

The equations of motion represent the first twelve elements of the system 
state vector. Hiese equations are summarized in Table 3.1 .1-1. 

The moment of inertia matrix I in general is given by 


I 

-I 

-I 

X 

xy 

zx 

-I 

I 

-I 

xy 

y 

yz 

-I 

-I 

I 

zx 

yz 
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for the moment axis terms, i.e., 1^, aind the product of inertia terms, 
i .e . , I 

zx 


Equations of 
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The aerodynamic force and moment coefficients and plume forces eure 
defined as functions of angle-of -attack, a, and angle-of -sideslip, 0, as 
shown in Figure 3. 1.1-1. The body referenced relative velocity vector, 
removing the wind velocity, v , from the vehicle velocity, is given by 


(B) B I (B) B LL (LL) 

v=v - Cv=v - C V 
-r — -w — -w 


(25) 


where v ^ ^^^ is the local -level referenced wind velocity vector. The 
following equations define a and 0 in terms of the components of v 


« = tan-^ (^) 
v 

^1 


(26) 


0 = sin ( ) 

V 

m 


(27) 


where 


V = (v ^ + V ^ + V ^) 
m r^ r^ r^ 


(28) 


The resultant thrust force f„ 

-T 


(B) 


-T 1=1 1 


T. 

1 

0 

0 


is expanded as 


~ i=1 ^i \ 


(29) 


B ^ 

where the transformation matrix transforms the magnitude of thrust 

for each thrusting device, f ^ , from its center-line to the body 

i 

coordinates. The general equation for f^ is 

i 


RELATIVE 

VELOCITY 

VECTOR 




-PA 

T. T. a. e 

1 ivac 1 


vacuum thrust 


P s atmospheric pressure at motor exit 
i 


A = motor exit cone area 
e 


B 

The matrix is different for the SSME’s and SRB' s and is given by 

f*B^MP MP^G- G^<L SSME 


C = transformation from the engine mount plane to the body 
coordinates 

MP G 

C = transformation from the gimbal reference plane to mount plane 
(structural deformation) 

G ^ 

C = transformation from enterline to the gimbal reference plane 
B 

C = transformation from SRB nozzle centerline to the body 


coordinates (gimbal angles ) . 




The resultant thrust torque is the sununation of the torque contribution 
from each thrusting device and is given by 


-T 


where 


iii 


, (B) 
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(B), B^<L 
r ) X C. 
-eg 1 


'T. 

1 

0 

0 


(31 ) 


(B) 

1 


body coordinates of the thrust reference point for the 
i^^ thrusting device . 


3.1.2 Measurement Equations 

Ihe measurements assumed available for the filter /smoother algorithm 
include inertial platform acceleration and attitudes, ground based radar 
tracking, SRB' s head pressure, SSME' s chamber pressures, liquid flow 
rates, pressurant flow rates. Ihe ET volumetric levels are available; 
however, due to their limited number (4), they may only be used for 
alternate checks of the filter /smoother algorithm perf ormamce . 

The propulsion related measurements will be treated in a separate 
section. In the following, the inertial platform acceleration measurements, 
attitude measurements and ground based tracking measurements models will 


be described for later linearization. 


3.1 .2.1 PlatfozTO Acceleration Measurements 


Accelerometers mounted orthogonally on an inertially stabilized 
platform, not located at the vehicle center of gravity, sense externally 
applied special forces and accelerations due to body rotation. The 
accelerometer measurement is modeled by 


(S) 
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-m 


S^P P^P’ 
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P' B 
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(B),, 
- r )J 
-eg 


+ b 
—a 


(S) 

+ V 
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(32) 


where 




p^p. 


P* 



= transformation from platform coordinates to sensing coordinates 

= transformation from misaligned platform coordinates to 
platform coordinates 

= transformation from body to misaligned platform coordinates 
= body coordinates of the platform center 
= accelerometer bias vector 


—a 

= accelerometer measurement noise vector 
-a 


3. 1.2. 2 Platform Attitude Measurements 

The inertially stabilized platform for the STS is a four axis IMU 
with a redundant roll axis [4 1. Vehicle body attitudes ere generated 
via quaternions [5]. It is assumed that an ecjuivalent representation 


can be made to obtain vehicle attitude by a three rotation sequence of 
roll, pitch, yaw to transform from inertial to body coordiantes . Ihis 
approach has been used in reference [6]. 

Ihe attitude angle measurement model is given by 


-m 


. (S) (S) 


(33) 


where 


P P’ 

^ = platform misalignment bias vector (used to formulate C ) 

(S ) 

y_ = attitude measurement noise vector . 

The transformation matrix used to treinsform from body to inertial 
coordinates in terms of the elements of the 6 vector is given by 


cos0cos4> sinq)sin6cos4» 

-cosq)sin(l< 


cos 9 sin 0 cos 4 » 

+sin<psin^ 




cos0sin<|> 


sinq>sin0sin4> 

+cos(pcos4< 


cosq>sin0sin4' 

-sin(pcosi|> 


-sin0 


sinipcos0 


cosq)cos0 


(34) 


3. 1.2. 3 Ground Based Tracking Measurements 

Ground based radar tracking devices can provide measurements of 
range, azimuth and elevation from the radar sight to the vehicle. Azimuth 
and elevation cure established relative to the sight's local level. If 
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the tracking device is a passive optical tracker (not laser) then only 
azimuth and elevation measiirements are available requiring more than one 
to establish position information. 

Defining Xi y, and z as the local easti north and up position of the 
vehicle relative to the ground based tracking device i the radar measurement 
equations &re given by 

ij 

p = (x^ + y^ + z^) + b + V (35) 

P P 

A = tan“^ (— ) + b + v (36) 

y A A 

E » tan ^ (■ — — ■ . ? ) + b_ + AE + v_ (37) 

.n — ^ 

Y X + y 

vdiere 


b^> b^, b^ a> range j azimuth, elevation biases 
AE = atmospheric refraction correction 


V , V,, v_ range, azimuth, elevation measurement noise, 
p A E 


Ihe position vector of the vehicle relative to the tracking device is 
given by 


A . (LL) LL^ECF ,ECF ECI (I) (ECF), 

‘ £ -£bdr ’ 


(38) 


z 


/ > 
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where 

LL GCF 

C > transformation from earth center fixed to local level 
ECF ECX 

C ■ eai -h centered inertial to earth centered fixed 
(ECF) 

r^j^ - position vector of tracking device in ECF coordinates. 

LL ECF 

The transformation matrix C is given by 





-sinX 

i 

-sinLcosX 

cosLcosX 


LL^ECF ^ 

■ cosX 

! 

-sinLsinX 

cosLsinX 

(39) 


! 0 

cosL 

j 

sinL j 



where L and \ aura the geodetic latitude and east longitude of the 

ECF GCJ 

device. The transformation matrix C is given by 


ECF^ECI 



sin(u) (t 
E 

■ ^RNP 


cos [w_(t 
E 

^RNP 

0 

0 



[RNP] (40) 


where 


u) =x earth rotation rate 

= ti"*® tag for RNP matrix 

KNF 


The position vector, 


(ECF) 


of the tracking device is given by 
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(ECF) 

^RDR 
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cos L 


+ (1 - e)' 


■ 2, 

sin L 


yj cos^L + (1 - e)^ sin^L 
R_(1 - e)^ 

( ' -" ■■■ 

\l 2 , ,, ,2 , 2 , 

j \ cos L + (1 - e) sin L 

I 


h) cosLcosX 


+ h) cosLsinX 


+ h) sinL 


where 

Rg = equatorial radius of Fisher ellipsoid 
e = flattening of Fisher ellipsoid 
h » altitude of the device above Fisher ellipsoid 


(41 ) 


I 



3 .2 Linearized System State and Measurement Eiquations 


The vehicle equations of motion are nonlinear fur tions of their 
motion variables and are implicit functions of other cements of the 
system states. Ibe measurement equations involve simil 2 ur function rela- 
tionships. The lineairizations for the filter/smoother algorithm require 

( B ) 

partial derivatives with respect to the motion variables, i.e., v emd 
6, and w‘th respect to other elements of the state vector, yielding 
explicit functional relationships for the elements of interest. 

For system state equations the partial derivatives will be presented 
in section 3.2.1 for the state elements in order of occurrence for the 
first twelve states. Other partial derivatives for ceuididate suate ele- 
ments will follow in section 3.3.1. The measurement equation partial 
derivatives for the first twelve states will be presented in section 3.2.2. 
Partial derivatives of the measurement equations for other candidate states 
will be pi esented in section 3 .3 .2 . 

Thti resulting partial derivatives are imbedded into the linearized 
system state matrix, F(x(t), t), as shown in Figure 3.2-1 . A corresponding 
linearized measurement matrix, H(^), is similarly formed with the 
measurement equations* partial derivatives. 

3.2.1 .System State Partial Derivatives 

Partial derivatives of each of the equations listed in Figure 3.2-1 
are developed in their order of occurrence with respect to the order of 











th« corrajponding states. Partial darivativaa of thrust tarms ara prasantad 
as though for a singla davloa. 


Inartial Position Rata Equation 

• (1 ) 

!2he first nonzezt) partial darivativa of tha r oquation is with 

^ ^ <B) 

rf^spact to V 






( 42 ) 


The second nonzero partial darivativa is with respect to This partial 

derivative results in a third order tensor and occurs frequently in later 
developments. The generalized form is presented in Appendix A. 


Body Velocity Rate Equation 

The partial derivative of with respect to r^^^ is given for 

altitude terms approximately as 




( 43 ) 


where 


8v Av - pAv 
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( 44 ) 
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The partial derivatives of 
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occur frequently and are given in Appendix B. 

The gravity vector partied derivative with respect 
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where 


H = gravitational constant. 

(B) 


0v 


The partial derivative, 


8r 


(I) 


, is the sum of the matricies in equations 


43 and 45. 


* ( B ) ( B ) 

The partial derivative of v with respect to v is given by 
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v^ere 
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{W')' - skew symmatric matxix made from the elements of the vector u> 
and equivalent to the cross product operator u x ( ) . 

• (B) 

The partial derivative of v with respect to 0 is 


• (B) 2 

5v pAv dv 0v pAv 0c^ _ 0v 

— m ^ m -r ^ m -f 0 a -r 
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The partial derivatives of are given in Appendix B and the partial 

-r 

0V 

•IT 

derivative is given in ^^pendix A. The last partial derivative 


is given in Appendix C. 
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The partial derivative of v with respect to (o is 


a*(B) , 2 - 

8 v pAv^ 3c^ 

a» "2111 ST * 


2 m 


(48) 



Euler Angle Rate Equation 


Die Euler angle rate equation is a function of both the Euler angles 
and the inertial ratec. The linearization will yield the two associated 
matricies . 

First with respect to the vector 0 , the following matrix results 


30 

30 


2 2 

qco89tcm0- rsin9tan0 qsinfsec 0 + rcos9sec 0 0 
-qsinf - rcos9 0 0 
qcos9sec0 - rsin9cos0 qsin9sec0tan0 + rcos9sec0tan0 0 


( 49 ) 


The partial derivative of 


0 with respect to to is 
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Inertial Angular Acceleration Equation 

The first partial derivative of this equation is with respect to the 
vector Using the approximation indicated in equation 43 , this 


partial derivative is 
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Next, with respect to the vector v , the partial derivative is 
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Hie partial derivative with respect to the vector 6 is 
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The final partial derivative for the first twelve states is with respect 
to the vector o). This operation yields 
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3.2.2 Measurement Partial Derivatives 


'Rie measurements assumed to be available > as discussed earlier i include 
grovuid based radar tracking i inertially stabilized platform attitudes 
relative to the vehicle body, and stabilized platform mounted 3 axis 
orthogonal accelerations. As with the state dynamics matrix, the measure- 
ment equations are linearized about the best state estimates. 


Radar Track Measurement Equation 

Referring to the radar track measurement equations, the required 
partial derivatives are 
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The rest of the required partial derivatives cure 
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Inertially Stabilized Platform Attitude Equation 

Ihe inertial platform is assumed to provide attitude angle measure- 
ments of the true attitude plus ein attitude bias plus measurment noise. 

The peurtial derivative of the measured attitudes with respect to the 
vector ^ yields an identity matrix. 

Accelerometer Measurement Equation 

The accelerometer senses specific body forces excluding gravity along 
the sensing axes. With reference to the accelerometer equation, the 
partial derivative with respect to r^^^ is 
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•Hie partial derivative with respect to v yields 
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For the partial derivative of the accelerometer wit! respect to 
the vector 0, the measurement equation is temporeirily rewritten as 
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where the vector represents the sum of the aerodynamic, thrust, 

P* B 

plume and rotational coupling terms. Ihe matrix C is the same 
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matrix as C . The required partial derivative results from 
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Ihe partial derivative on the right hand side is developed in Appendix A 
(B ) 

with the vector s representing the su.n of the terms indicated above . 

The final partial derivative for the accelerometer measiirement is 
with respect to the body rotation vector u. Defining 
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3.2.3 Additional Parameter Partial Derivatives 


The mathematical developments are presented in this section for the 

partial derivatives of the system and measurement equations to allow for 

additional candidate parameters to be included in the estimation algorithms. 

These parameters include center-of -gravity, r , moments of inertia, I, 

”^g 

wind velocity, v^, and inertial platform tilt errors. Aerodynamic and 
plume parameter partial derivatives aure also presented. 

The computer program is being structured to permit these parameters 
to be easily incorporated without significant impact on the program code. 


3. 2. 3.1 Center-of -Gravi ty 

From equation 23, the partial derivative of angular acceleration with 

respect to r is 
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From equation 32, the partial derivative of the measured acceleration with 

lespect to r is 
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3. 2. 3. 2 Momenta of Inertia 

The momenta-of-inertia are grouped into "principal" terma, and 

croaa product terma > i^p* From equation 24, theae vectora are defined as 
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With these definitions, equation 23 is rewritten as 
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\vhere IT repi ;ents the sum of the nonrotational torques in equation 23. 
Defining an intermediate vector a as 
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the partial derivatives of the angular acceleration with respect to 
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3. 2. 3. 3 Wind Velocity 

From equations 20 and 25, the partial derivative of the vehicle 
acceleration with respect to v is 
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The first of the pfurtial derivatives in equation 79 can be obtained 
from the following equation 
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From equation 80, the following is obtained 
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Denoting the elements of the matrix C as ^ , c^ 2 » ®tc., the following 
equations eu:e obtained for the partial derivatives of o and B with 
respect to v ; 
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3. 2. 3. 4 Inertial Platform Tilt 


Temporarily rewriting equation (32) as 
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where 


66 = vector vghose elements are the axes misalignments 

p* B 

s = sura of the bracketed terms in equation 32 multiplied by C . 


Hie following peurtial derivative of the measured acceleration with 
respect to 66 is obtained 
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3. 2. 3. 5 Aerodynamic and Plume Parameters 

A linear model for the aerodynamic and plume characteristics is 
used. Ibis model is e^qaanded as 
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where additional terms to represent rates, cross couplings, and controls 
can be included. 

The laasic approach of establishing the partial derivatives will be 

illustrated for a couple of terms, c, and c • Using these example 

— la ~ma 


illustrations , the rest of the candidate parameters can be similarly 
obtained. From equation 20, the following peurtial derivative is obtained 
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vrtiere 


[U] = unit 3x3 matrix with one's (1 ) on the diagonal and zeros 
off the diagonal 

From equation 23, the partial derivative of angular acceleration with 

respect to c is 
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The corresponding partial derivative with respect to c^ is 
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Ihe static aerodynamic coefficient model has been obtained by a multiple 
regression analysis of the current aerodyneunic tabular data. Ihis model 
is presented in Appendix D with the associated regression coefficients. 



3 .3 Propulsion Parameter States and Measurements 

A candidate approach for incorporating the NASA propulsion model's 
capabilities has been identified. This approach utilizes nominal pre- 
dicted values of thrust > pressure, propellant and pressurant mass flow 
rates, and utilizes sensitivities or partial derivatives of these 
variables with respect to the independent parameters selected for esti- 
mation by the algorithm. 

The approach is to include deviations f • om nominal values of 
measured chamber pressure, power level, propellant and pressurant mass 
flow rates as states. Thj models as?; tad for these deviations are time 
correlated random processes. Then as states, partial derivatives of 
the first twelve states with 'espect to these variables will be required. 

For the SSME and SRB, this modeling approach is discussed in the 
following. Additionally, the necessary partial derivatives of the first 
twelve state variables with recepct to the additional states are presented. 


3.3.1 SSME Propulsion Parameter Model 

For the SSME, the total actual values of vacuum thrust and oxidizer 
mass flow rates are modeled by 
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Ihe measurements of fuel mass flow rate, pressurant mass flow rates and 


power level are modeled as 
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These measured quantities include measurement noise s ^ ^ and potential 
bias states b^ ^ modeled as random constants. In these measurements) 
the A'd variables are to be included as states in the estimation 
algorithm. If the nominal values are zero or unknown, then the A'd 
variables absorb the entire estimate. Where required, the estimate for 
the variables used in the estimation algorithm is formed using the 
nominal and the estimate of the deviation, etc. In example, thrust and 
fuel mass flow rate estimates are formed as 
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The deviation or A'd measurement variables are modeled as time 
correlated random vairiables. *niis permits these variables to vary within 
a band of frequencies . The typical model is then given as 

ir *(> = -- A()t- s,, (100) 

dt t, , ^( ) < > 

vrtiere the parenthesis ( ) would be replaced by the variables, i.e., m . 

* 2 

For the SSME, an additional variable Ac is modeled as in equation 

mult 

100 and included as a state variable with the A'd measurement 
variables . 

Hie thrust deviation is expeinded as in the following truncated Taylor 
series as a fianction of the independent parameters. 
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In the V and w equations, with equation 101 replacing f^ 
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the partial derivatives of f^ with respect to the A'd variables are 
obtained directly from equation 101. 



It Is desirable to include vehicle mass bias as a state. The SSME* s 


system contribution to the mass deviation is given by 
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In equation 101, the Am contribution to the mass deviation is not 

2 

available from measurements. As with the thrust deviation, this quantity 
is formed as 
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which is in terms of other estimated state Vcuriables. In equations 101 
and 103 the deviation in mixture ratio, AMR, is obtained algebraically 
from 
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The partial derivatives for the SSME above have been incorporated into 
the estimation algorithm as functions of engine power level . 


3.3.2 SRB Propulsion Parameter Model 


Ihe approach for the SRB modeling follows closely that used for the 

SSME. Candidate independent parameters include propellant burn rate 

exponent, a, and motor efficiency coefficient, c . Others can be added 

m 

using this technique. 

The actual value of vacuum thrust is given by equation 92 . Ihe 
only measurement available for the SRB is the total pressure at the 
forward head end of the motor case and is modeled as 
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vrtiere b^ ^ and s^ ^ represents a bias auid me-3urement noise respectively. 

The independent parameters, Aa and Ac , are included in the 

ret 

model as states. The model assumed can be as given by equation 100 or 
another suitable dynamical process, i.e., random constant. 

The thrust deviation is given by the following truncated Taylor 
series as a function of the cauididate independent parameters . 
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The partial derivatives for the and w equations with respect to 
the independent parameters are obtained directly from equation 106 . The 
mass deviation equation for the SRB is given as 
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where 
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X oa 

The head pressure deviation, AP , is expanded similarly 

APo = ~~ Aa + . . . (109) 


A simplified model for the SRB's thrust, head pressure and mass flow 
rate has been developed that can be directly incorporated within the filter 
algorithm for estimating burn rate coefficient, nozzle coefficient, mass 
flow rate, etc. This model, to be described below, uses apriori specified 
burn area and port volume as a function of burn depth into the propellant 
grain. Prom this simplified model analytical partial derivatives required 
by the estimation algorithm can be obtained. 

The thrust is given by 
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c = nozzle coefficient 
rr. 

c^ = thrust coefficient 
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c = characteristic exhaust velocity 
m = mass flow rate 

TWo of the required partial derivatives with respect to mass flow rate 
and nozzle coefficient are easily obtained, vis 
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Ihe partial derivative with respect to burn rate coefficient is 
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where it has been assumed that c is not a function of a. Using the 
"ideal" expression for c^ [7] 


f Y>1 T+1 


1+1 


-1 


[1 - 


(_e) 

P_ 




where 

Y = ratio of specific heats 

= motor nozzle exit pressure 

» ambient atmospheric pressure at nozzle exit 
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the first partial derivative in equation 113 is 
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Ito evaluate the second partial derivative in equation 113, the following 

equation for pressure [8] is used: 
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where 


p » propellant density 
P 

= propellant burn area 
n = propellant burn rate exponent 


Ihe following partial derivative is then obtained 
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The last partial derivative in equation 113 is obtained from 
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The resulting partial derivative is 
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Tb utilize the head pressure measurement and its sensitivity to 
parameter variations, the following equation [7] is used 
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where 

R > gas constant 
T a gas absolute temperature 
c • port circumference 
Ap » port cross section eurea 
- propellant burn rate 

t 3 distance from motor nozzle to pressure measurement point 


This equation assumes a cylindrical c with an approximately constant 
cross sectional area. 

The partial derivatives with respect to burn rate coefficient and mass 
flow rate are 
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In aquations 121 and 122 a cylindrical port has b«an assumed in determining 

the port volumn V^. Equation 122 was obtained from equation 120 by 

replacing the term c r. p 1 by m. Ihe partial derivative of P. 

op 

with respect to c is obviously zero. In using these analytical 
m 

partial derivatives ^ the basic performance measures of tiirust; mass flow 

rate I head and nozzle pressures, etc. are matched between this model and 

the NASA SOBER internal ballistics routine results. The burn area and 

port volume are adjusted in the simple model to obtain the agreement. 

Ihen using the adjusted area and volume as a function of burn depth, the 

partial derivatives are evalauted. 

Ihe inclusion of solid propellant thickness, x, as a state variable 

necessitates the development of the partial derivatives of x with respect 

to the solid propulsion parameters and the partial derivative of the 

measurement P with respect to x. These partial derivatives are given 
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as 
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Here I the partial derivatives of c , A^, and A^ with respect to x are 
evaluated numerically. 

Finally, the partiad derivative of thrust with respect to x is 
given by 
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3.3,3 Vehicle Miss 


origi:ial 

OF POOR 


V- 


The total rate of change of vehicle mass is given by 
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The first tvro terms in this equation are the apriori assumed nominal 

values. The third and fourth terms were discussed earlier. Ihe last 

term should be zero; howover it can include a mass bias uncertainty 

The equations, state and measurment, in which mass occurs are the 
( B ) 

and ^ equations. Assuming equation 123 can be summarized as 

• • 

m + then the mass can be written as m + Replacing this 

expression for the mass in the two indicated equations yields the 
following peirtial derivatives with respect to the 
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4.0 PROJECTED ACTIVITIES DURING UPCOMING MONTO 



During the upcoming month the computer programs will be exercised 
on the NASA computers , and documentation of the computer routines will 
be developed. 

Due to minor delays, SDI is requesting a one month extension at no 
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APPENDIX A 


PARTIAL DERIVATIVE OF THE VECTOR 




V wrt 


0 


This partial derivative is one of several that occurs frequently in 
the formulation of the linearized system state and measuremetn eqautions. 
The desired peurtial derivative is 
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1 -k:os<pcos 4' V -sinipcos(|« 3 


(-sin0)v^ + (sin9cos0)v2 + (cos9cos0 ) v^ 
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The resulting matrix is given in Table A-1 . 





i 


APPENDIX B 


0v 


m da 


(B)’ 


dv 


m 


av av' av 


as aa aa , ... „ 

W riST- 3h- ah "ah 


These partial derivatives occur frequently and will be developed 
in this appendix. The equation for v is ' 


(B) B LL (LL) 
V = V - C V 
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Since the wind velocity i is only a function of altitude then 
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Ihe equation for — r— is 
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The following equations define the last three 


required partial derivatives 
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PARTIAL DERIVATIVE OF THE VECTOR 


V wrt 8 


The third of the frequently occurring required partial derivatives 


is 


cos6cos(t( V 
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The resulting matrix is given in Table C-1 . 
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APPENDIX D 


AERODYNAMIC MODELING REGRESSION ANALYSIS AND RESULTS 


Ihe aerodynamic data tables provided as IVBC3 data has been incorporated 
into an aerodynamic coefficient polynomial model. This modeling effort 
reduces the dimensionality of the munerical tables to one and reduces the 
storage requiranents for the aerodynamic model. 

The coefficient model used for the two stages differ slightly as a 
result of the available data. Ihe regression analysis led in the selection 
of the form of the aerodynamic model. Terms with insignificant correlation 
were eliminated from the model. 

In equation form, the first stage static coefficients of axial force, 

C ; normal force, C ; pitching moment, C ; rolling moment, C. ; side force, 

C^; and yawing moment, C^; are given below 




^ ^ ^A 2 ^ ^A 7 

a a ap p 


D-1 




+ C. 


* „2 

o3 


D-2 


c 

m 


C 

m 

o 


+ C o + C 


m 






D-3 



+ C- B + C otB + C. a 6 
B aB a B 


D-4 


ORIGINAL FAC 'I uii 
OF POOR QUALlTV 


- C„ + C„ 0 + C„ oft + C„ . a 0 


^ Yo 'Yg 


O0 


V 0 


D-5 


C«C +C 0+C a0+C« o 0 

n n n- n * n 2. 

o 0 a0 a 0 


D-6 


Ihe corresponding second stage model is given by the following equations 
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For the first stage, data from an angle-of -attack range of -6 to +6 degrees 
was used in the regression analysis. Data from a range of -8 to +4 degrees 
was used for the sec' id stage. Ihe results, TCXX,.., from the regression 

analysis is presented below for each of t!ie coefficients, C , above. 
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TABLE I>-2. First Stage Normal Et>rce Coefficient 
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TABLE D-7. Second Stage longitudinal Coefficients 


OF POOR QUALi 


v-i 

u 


M M n « 

0 c o o o 

1 I I I I 

UJ 

a O' 

C'i f5 O fN O 

c ^ IS cc >0 r< 

ri ^“5 o M sj 'C 

<I i-( Hi O'! IH r-i 

O 35 ^ O rr M 

I- »«i ->0 tv' O «■ 



n 

M 

<r 

«r 11 


■<T Ki 



0 

1 

C 

i 

0 

1 

c o 


O C 

0 

1 


UJ 

U! 

1 

UJ 

U: UJ 


1 1 
UJ UJ 

1 

UJ 


^3 



O' Tv 


O' O' 

r3 

C-J 

03 

rs 

'C 

C 4 « 


IS -<3 

0? 


tr 

hO 

03 

O L3 

(M 

X— 4 IS 

<■ 

rs 

o 

Os 

O 

CC 

<C 

IS <r 

O' 

z 

r*5 

L~ 

O' 

t*i <! 


11 c 

o 

u 

t*3 

«r 

sD 

r> IS 

s~ 

03 O' 

r*5 

H- 


cc 

33 

r-j O' 

U 

fs 

'C 


ti 


f3 

M C-i 

►- 

▼S .-i 

00 


c o o o c 


o o c o o 

I I 


0 o o 

1 



rs c-i 

c-i 

c-i rs 









rj 

rs 


O c 

o 

c o 


C 

o 

o 

c 

o 


e 

c 

o 


UJ 

Ui 

UJ 

UJ UJ 


UJ 

Lu 

Ll 

UJ 

UJ 


LU 

LU 

UJ 


c 

S3 

O' 

o c 

<E 

f-i 

rt 

c 

cc 

K! 

<c 

IS 

r3 


<i 

IS 

tT 

r' 

T r<3 

c-i 


CO 

c 

r-4 

*r 

O'! 

04 


t3 

CJ 

13 

03 

C' 

TS -9 

2 

•c 

C-4 

c 



r 

•C 

IS 

X 



rs 

'T 

r*v ^ 

u 

O' 

<♦ 

t3 

S3 

Tv 

u 


11 

r-4 

a 

fS 

S3 


to C' 

1- 

IS 


rs 

•<r 

r. 

h" 

f-J 

m 

O' 

H- 

c 

r-j 

rs 

r*3- 33 


r*3 


03 

O' 

r-4 


X 

IS 



03 

Tv 

rr 

S3 IS 


IS 

•r 


c 

c 


o 

S3 

05 


W 

<T 


rr <■ 


Cn 

O'! 


r-j 

C-i 



X 

S3 


o 

o 

c 

0 

0 


c 

O 

o 

o 

O 


o 

o 

O 


I I I I I III 


IS f-i «C O' 

h; T'. N •*-( UT 

O 111 «■ *0 i> 

r>l 'C 03 C Tx c 

» S 3 03 r -4 >0 

n a IS 03 53 r*3 t' 

U. I— IS S 3 >C IS S 3 

«X »-i »s rs ts 

CJ . . . . . 

O O c o c 


^ r *3 rs c-i ii3 

■c L"; rs o. f j 

O 0^ ^ t*3 

Cm' r-i C C xC o 

z o; CO CK IS o 

r -4 O f s -C •«? M O' 

U. t— CC IS S3 S3 >0 

z o o o o c 

U 

c o o o o 

I I I I I 


S3 O tS 

33 ^ O' 

C'i tn 

O <• M "C 

C4 ^ rj fs 

w r 1:3 IS rs 

U. U ';3 ^ 

Z I- o o o 

o . . » 


o o o o c 


I- I- K 

Ul UJ u 

CI3 c<3 cn 


z 

UJ 


z 

UJ 4^ 


*X O iMX O uliX O 

CUUTnOO* CUU3L3S50- OUUTpOO* 

^2 • ♦ * • ^3 xX ' ♦ • • ^3 ♦ ^3 

Ci.ZM'<TS3C0^ u.ZM<rs3 03rs CLZr3VSiQ3»'< 


03900000 -0.5A07i4?lE-02 0 . 8928n677E >04 

03719047 -0.50357133^-02 0 . 1 1 30953?.f.-03 


TABLE D-8. Second Stage Lateral Coefficients 
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